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Abstract 

The inclusive decay B — > X u £u is of much interest because of its potential to constrain 
the CKM element \ V u b\. Experimental cuts required to suppress charm background restrict 
measurements of this decay to the shape-function region, where the hadronic final state 
carries a large energy but only a moderate invariant mass. In this kinematic region, 
the differential decay distributions satisfy a factorization formula of the form H ■ J <S> S, 
where S is the non-perturbative shape function, and the object H ■ J is a perturbatively 
calculable hard-scattering kernel. In this paper we present the calculation of the hard 
function H at next-to-next-to-leading order (NNLO) in perturbation theory. Combined 
with the known NNLO result for the jet function J, this completes the perturbative part 
of the NNLO calculation for this process. 



1 Introduction 



The inclusive decay B — > X u iv is of much interest because of its potential to constrain the 
CKM element \V u b\. Due to experimental cuts required to suppress charm background, 
measurements of this decay are available only in the shape-function region, where the 
hadronic final state is collimated into a single jet carrying a large energy on the order of 
m&, and a moderate invariant mass squared on the order of m;,AQCD- Much theoretical ef- 
fort has been put into establishing a factorization formalism which enables the calculation 
of differential decay rates in this kinematic region. Early work in QCD was based on dia- 
grammatic approaches [1,2], whereas more recent papers [3-5] are based on soft-collinear 
effective theory (SCET) [6-8]. The main result of these works can be summarized in the 
following factorization formula for an arbitrary differential decay rate: 

dT~H-J®S, (1) 

where the symbol (g> denotes a convolution. The perturbative information is contained 
in the hard function H, which is related to physics at the hard scale m^, and the jet 
function J, which is related to physics at the intermediate scale m fe AQ CD . The object 
S is a non-perturbative shape function describing the internal soft dynamics of the B 
meson [9,10]. The factorization formula is valid up to corrections in which 
have been studied in detail in [11-13]. The hard and jet functions to next-to-leading 
order (NLO) in perturbation theory have been known for some time [3,4], and the jet 
function at next-to- next-to-leading order (NNLO) was obtained in [14]. 

The main purpose of this paper is to complete the perturbative part of the NNLO 
corrections to the factorization formula (CD) by obtaining the hard function to this order. 
The organization is as follows. In Section [2j we briefly outline how to obtain the hard 
function through a matching calculation in SCET. The task is to extract three Wilson 
coefficients Cj, which arise from integrating out the hard scale m& by matching the semi- 
leptonic b — > u transition current from QCD onto SCET. The discussion there makes 
clear that the principle technical challenge is to calculate the two-loop QCD corrections 
to the b —>■ u current. This loop calculation is the subject of Section [3l where we ex- 
plain our calculational procedure and give explicit results in terms of a set of harmonic 
polylogarithms. The method relies on a reduction to master integrals through integration- 
by-parts relations, which are then solved using differential equations. In Section HI we use 
our results to obtain the Wilson coefficients at NNLO; a phenomenological analysis of 
partial decay rates and the impact on the determination of \V u b\ is in progress and will 
be presented in future work. We conclude in Section [51 

2 The hard function in SCET 

The QCD effects in inclusive semi-leptonic B decays are contained in the hadronic tensor 
W^, from which any differential decay distribution can be derived. It is defined as the 
discontinuity of the forward matrix element of the current correlator T 7 ^, which is the 
time-ordered product of two semi-leptonic b — > u currents, J M = «7 M (1 — 75)6: 

W" v = ^Im ( B ( V ^M , T» v = i J d 4 xe iq - xr T{J^{S)r{x)} . (2) 
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Here q is the momentum carried by the lepton pair and v is the velocity of the B meson. 
Using the SCET formalism it is possible to show that the hadronic tensor obeys the 
factorization formula 



i,j=l ^ ' 



(3) 



We have introduced the vector p = m b v — q, which in the parton model is the momentum 
of the final-state jet into which the b quark decays, as well as its light-cone decomposition, 

= (n ■ p) — + p\ + (n ■ p) — = p% + pt + p\ , (4) 

where n and n are two light-like vectors satisfying n ■ n — 2. The object Hij is defined as 

Hij(n ■ p) = Ci(n ■ p)Cj(n ■ p) , (5) 

where the Wilson coefficients C{ arise from matching the semi-leptonic b — > u current 
from QCD onto SCET. In position space and to leading order in the heavy-quark limit, 
this matching is of the form 

3 

e- imbV - x u{x)^{l - 75 )&(a0 = T I ds Ci{s)x{x + sn)T^n{x^), (6) 



2^ / ds d(s)x(x + sn)TfH(x-), 
i=i J 



where we have followed the SCET conventions of [4]. The are a set of three Dirac 
structures, which we shall choose as 



n 



ft 



r? = 7 ^(i- 75 ), r^ = ^(i + 75 ), n = — (i + 75 ). (7) 

In practice, the matching calculation is carried out in momentum space and yields results 
for the Fourier-transformed coefficients, which read 

Ci(fi-p) = J dse lsn - p C t (s). (8) 

The matching coefficients are obtained by evaluating UV-renormalized matrix elements 
of both sides of ([6]), corresponding to calculations in full QCD and SCET. The calculation 
is simplest when the external states are chosen as on-shell quarks and both UV and IR 
divergences are regulated in dimensional regularization in d = 4 — 2e dimensions. In that 
case the loop corrections to the SCET matrix elements are given by scaleless integrals and 
vanish, so that the result is just its tree-level value multiplied by renormalization factors 
from operator and wave-function renormalization. The QCD result is written in terms of 
three Dirac structures multiplied by scalar form factors, which we shall define according 
to 

(u(p)\J»\b(p b )) = D l u(p)Y(l-l 5 )u(p b ) + D 2 u(p)^(l + l5 )u( Pb ) 

m b 

+ D 3 w(p) — (1 + 75MP6), (9) 



2 



where u(p) and u(pb) are on-shell spinor wave functions, pb and p are the momenta of the 
b and u quarks respectively, and p 2 = 0, pi = m b . We shall always work in the reference 
frame where the perpendicular components of the external momenta vanish, and where 
p^ = m b v^ and p M = (n • p)n tl /2. Then the three Dirac structures multiplying the Di 
correspond to those in ((7|) in an obvious way 

To determine the Wilson coefficients Ci we also need the SCET matrix element, for 
which we can make an important simplification. In general, the result involves a renor- 
malization matrix Zij applied to the bare SCET current operators. However, we can use 
that the partonic expression for the quantity J ® S in the factorization formula ([3]) for the 
hadronic tensor is independent of the coefficients Hij that multiply it. This implies that 
the operator renormalization matrix is just the unit matrix multiplied by a single scalar 
factor Zj. Moreover, for on-shell matching the wave function renormalization factors 
in SCET are unity, and the SCET spinor wave functions correspond to those in QCD. 
Therefore, the coefficients C, can be obtained through the relations 

Ci(n-p) = lim Zj 1 (e,m b ,n-p,{i)D i (e,m b ,n-p,n) (2 = 1,2), 

C 3 (n-p) = lim Zj 1 (e,m b ,n ■ p, /i)-^-^-D 3 (e,m b ,n • p, fi) ■ (10) 

The renormalization factor Zj can be determined in two different ways. The first is to 
require that the matching relation ( jTOl) is free of IR poles in dimensional regularization, 
which allows one to deduce the UV structure of the SCET currents from the IR structure 
of the Di. A second method is to determine the UV poles of the object J <g> S in the 
parton model, using the two-loop anomalous dimensions for the jet and soft functions, 
calculated in [14] and [15,16]. Agreement between the two methods is an important check 
on the factorization formalism, and also on the two-loop calculation of each function. The 
agreement will be verified in Section H] below. 

We end this section by pointing out a subtlety in the matching calculation related to 
heavy-quark loops, which first becomes relevant at NNLO. Whereas the partonic matrix 
elements in QCD are calculated as an expansion in a s in the MS renormalization scheme 
in a five-flavor theory, where n/ = n\ + nn with nu — 1 for the b quark, in SCET 6-quark 
loops are absent and the matrix elements are calculated as an expansion in a four-flavor 
theory. To match results in the two theories as in ( ITUl) . it is necessary to express the UV 
renormalized results in five-flavor QCD in terms of the four-flavor parameters of SCET. 
To achieve this, one renormalizes the coupling constant in the rif = nh + ti\ flavor theory 
according to a^ are = Z™ h+ni a si with (see e.g. [17]) 



>yn h +ni _ 1 _ 



11 4 4 

—C A - -T R n f + -T R n h (l - N e 



'111 



The function N e is fixed such that a s is the MS-renormalized coupling in the four flavor 
theory. Its value is 

N{e)=eJe {^J T{l + e) - (12) 
Results for the scalar amplitudes Di in this renormalization scheme can be obtained from 
those in the MS scheme in five-flavor QCD by making the replacement 



a s — > a s [ 1 + ^3 T R n h 



L + e\ L -\ + e + —L 

2A V 3 12 6 



+ ..., (13) 
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where L = lia.fi/mi,. After applying this decoupling to the Di, dependence on rih in the 
pole terms, and thus Zj, drops out. This must be the case, since heavy quark loops do 
not exist in SCET, where the h quark field is treated as in HQET. This same procedure 
was used in the completely analogous case of matching the b — > s current at q 2 = in [18]. 

From the above discussion, it is obvious that the main technical obstacle to obtaining 
the Wilson coefficients Cj is the calculation of the QCD form factors Di. This will be the 
subject of the next section. 

3 Two-loop QCD corrections to the b — > u current 

In this section we perform the calculation of the renormalized scalar form factors Di at 
two-loop order. We begin by outlining the calculational procedure in Section 13.11 and 
then give the final results in Section 13.21 

3.1 Calculational procedure 

In this section we describe some technical details involved in obtaining the two-loop QCD 
corrections to the b u current. The main task is to evaluate the bare two- loop amplitude 
by calculating the Feynman diagrams in Figure [H This bare amplitude contains both UV 
and IR divergences. The UV divergences are removed by counterterms related to b and it- 
quark wave-function renormalization (on-shell scheme), coupling constant renormalization 
(MS scheme), and mass renormalization (on-shell scheme). 

The calculation of the individual two-loop Feynman diagrams proceeds as follows. 
First, by doing tensor decomposition, we extract the contributions of each diagram to the 
form factors in 0. At this level, these contributions are written as linear combinations 
of certain scalar integrals. Second, this rather large set of scalar integrals is reduced to a 
much smaller set of master integrals using the Laporta algorithm [19], which is based on 
the integration-by-parts identities introduced in [20,21]. A very useful tool for performing 
this reduction is the integral reduction program AIR [22], written in Maple, and we have 
used this program in our calculation. 

A typical master integral depends on rrib, the dimensionless variable s = (pb—p) 2 /™^ 
and the parameter e = (4 — d)/2 of dimensional regularization. Some of the simpler master 
integrals (those with three or less propagators), are easily solved using the standard 
technique of Feynman parameterization. In most cases, it is straightforward to obtain 
exact results in e, which involve hypergeometric functions or their generalizations. These 
can be expanded around e — > using the Mathematica program HypExp [23,24]. For the 
more difficult master integrals, we have used the differential equation technique [25] (for 
a recent review, see [26]). This involves solving a set of differential equations obtained by 
differentiating the master integrals with respect to the variable s. The solutions to the 
differential equations determine the master integrals as a Laurent series in e, up to their 
values at the boundary point s — 0. In some cases, these constants can be determined 
by requiring that the coefficients in the Laurent expansion are finite in the limit s — > 0. 
In other cases, there is no choice but to calculate the e-expansion of the two-loop master 
integral at the point s = 0. The solutions to the differential equations involve the harmonic 
polylogarithms (HPLs) introduced in [27]. For their numerical implementation and also 
some symbolic manipulations, we used the Mathematica package HPL [28]. 
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Figure 1: Two- loop corrections to the b — *> u left-handed current. The incoming 6-quark 
and the outgoing w-quark are represented by thick and thin solid lines, respectively, while 
dashed lines represent gluons. Fermionic bubbles with 6-quarks and lighter quarks (the 
latter being treated as massless) are shown by thick and thin circles. Diagrams where the 
light fermionic bubbles are replaced by gluons and ghost-particles are not shown explicitly, 
but they are taken into account. 



We have checked our results in several ways. First, we have used the numerical method 
of sector decomposition [29] to evaluate the master integrals for various values of s, and 
checked that they agree with the analytic results. For this we have used self-written 
code, and also the publicly available C++ program described in [30]. Second, we have 
obtained results as a double series in e — > 0, s — > using two different techniques. One 
is to expand each master integral as a series in s —>■ before doing the loop integrals 
using sector decomposition, the other is to obtain results for each diagram at s = and 
then recover the s-dependence using differential equations. We then checked that these 
agree numerically with the expansion of the analytic results in the same limit, up to the 
first five or six terms around s — > 0. Finally, we were able to transform our basis of 
master integrals into that used for the two-loop calculation of the vertex corrections in 
B — > 7T7T, presented in [31,32]. For some of the master integrals, we used these results to 
help convert numerical results for the boundary conditions into analytic results in terms 
of constants like tt. 

To illustrate the method of differential equations in our application, we take as an 
example the first diagram in the second row in Figure [U In this case we have four master 
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integrals hi, h 2 , h 3 and h 4 , reading 

, d d i d d r 1 

hi(s) 
h 2 (s) 



(2ir) d {27i) d [{£ + p b f - m 2 ] [(/ + r + pf] \{r + pf] 
d d £ d d r 1 



(27r) d {27c) d [{£ + p b f - m 2 ] [(£ + pf] [(Z + r + p) 2 ] [(r + p) 2 ] ' 
d d £ cZ d r 1 
(2^p (2^p p + p 6 ) 2 - m 2 ] [r 2 ] [(/ + r + p) 2 ] ' 

n /" _d d ^^ d r_ 1 

4lSj 7 (27r) d (2vr) d [(£ + ^) 2 - m 2 ][r 2 ][(£ + p) 2 ][(/ + r+j9) 2 ] ' 1 j 

They satisfy the differential equations 

dhAs) 



ds 







dh 2 (§) 1 (4cZ + Ads - 16s - 12) , _ 1 (-3d + 8) , _ 

— ; = ; ; hoys) 7, — : thus), 

ds 4 s(l -s) V ; 4 m 2 s(l - s) V ; ' 

^Mf) 1 (4-*) u f * 1 m 2 (l-g)(d-4) 

tZ/i 4 (s) 1 (3d + 5ds-8-20s) 1 (-3d + 8) 

" ~4 1(1=1) Ms) "4mgS(l-S) h ' 3is) - (15) 

Obviously, fti has to be calculated using the standard technique of Feynman parameteriza- 
tion. The dependence of h 2 on s can then be determined by solving the second differential 
equation, in which hi plays the role of a given inhomogeneity. The requirement that h 2 
is non-singular for s — > uniquely determines the function h 2 (s). The s dependence of 
the functions ha and hi can be obtained by solving the corresponding two differential 
equations simultaneously (as an expansion in e). Specifying h 3 {s = 0) by means of stan- 
dard Feynman parameterization and imposing the additional requirement that h±(s) is 
non-singular for s — > 0, uniquely determines h 3 (s) and h^s). 



3.2 Renormalized scalar form factors 



We now give results for the UV-renormalized form factors in (j9]), which we expand in ct s 
according to 

We start by listing the results of the one-loop contributions. To this end we further 
decompose the quantities Df' as 

nd) 



d(1) d(1) 



(16) 



The Laurent expansion coefficients of the poles and constant term have been known for 
some time [7], whereas the terms proportional to e and e 2 are new. Note that terms up 
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to e 2 are needed to correctly extract the Wilson coefficients Ci through (JTU1) . The explicit 
results for the I&> in (fT6l) read (recall s = (pb — p) 2 /ml) 



= -1 



L (-2),l 
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^(-i),i ~~ ---2L-2F 4 



F 4 



= -6-5L-2L 2 -^-3F 4 -4LF 4 + ^-2F 5 -4F 10 

= -12-12L-5L 2 -^-^-^-8F 4 -6LF 4 -4L 2 F 4 -^ + 
W' 1 3 24 6 6 

4_4 + 2L^ _ ^ _ + _s _ ^ _ g ^ _ + 2^0 _ _ 

s s s s 

C(3) 



4F 13 - 8F 17 + 
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10L 3 2L 4 7T 2 5Lvr 2 L 2 rr 2 tt 4 



i?S , = -24 - 24L - 12L 2 - — — - — 16F 4 - 

( 2 )' x 3 3 2 12 6 160 

o 8L 3 F 4 tt 2 F 4 1 „ 8F 4 8LF 4 2L 2 F 4 

16LF 4 - 6L 2 F 4 — i -± - -Ltt 2 F 4 + + — — + — — + 

3 4 3 s s s 

— - - 8F 5 - 6LF 5 - 4L 2 F 5 - — ± + + — ^ - 3F 6 - 4LF 6 + 
12s oss 

^-2^-16^-12^0-8^0-^ + ^ + ^-6^- 
s 3 s s 

2.F 1 2.F 1 
8LF n + — - - 4F 12 - 6F 13 - 8LF 13 + — — - 4F 14 - 4F 15 - 12F 17 - 

s s 

AF 5^(3") 2 2 

16LF 17 + — il - 8F 18 - 8F 19 - 8F 20 - 16F 21 + + -LC(3) + -F 4 C(3) 

s 6 3 3 

_ p(i) _ n 

rL (-2),2 ~~ ^(-1),2 ~~ U 

(!) 2 ^ 2F 4 2F 4 

(°)- 2 ~ s ~~ s 2 s 

(1) 4 4L _ 2F 4 _ 4LF 4 2F 4 4LF 4 2F 5 2F 5 4F 10 4F 10 

y >' s s s z s z s s s z s s z s 

(1) 8 8L 4L 2 tt 2 4F 4 4LF 4 4L 2 F 4 vr 2 F 4 4F 4 4LF 4 

( 2 ) > 2 s s s 6s s 2 s 2 s 2 6s 2 s s 

4L 2 F 4 | tt 2 F 4 2F 5 4LF 5 | 2F 5 | 4LF 5 2F 6 | 2F 6 4F 10 

s 6s s 2 s 2 s s s 2 s s 2 

SLFi n 4-Ti n SLFi n 4-Ti i 4.r i i AFi 3 4.r i 3 8-f i •? &Fi 7 



+ — - + — 1 - + — : ;=F + 



s 2 s s s 2 s s 2 s s 2 
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d(i) _ d(i) _ n 

rl (-2),3 ~~ (-1),3 ~~ U 

„(i) _2 2F\_AF\ 



(°)' 3 s s 2 s 

r{1) _4 _ 4L 2^ 4L^4 _ 10^4 _ 8L^4 | 2F 5 4F 5 | 4F 10 8F 10 

W' 3 s s s 2 s 2 s s s 2 s s 2 s 

(1) 8 8L 4L 2 vr 2 4F4 4LF 4 4L 2 F 4 tt 2 F 4 20F 4 20LF 4 

( 2 ) ' 3 s s § 6s s 2 s 2 s 2 6s 2 s s 

8L 2 F 4 tt 2 F 4 2F 5 4LF 5 10F 5 8LF 5 2F 6 4F 6 4F 10 

^ a 1 ^9 ^9 ^ ^ ^9 ^ ^9 

s 3s s 2 s z s s s z s s z 

8LF W 20F 10 1QLF 10 | 4F n 8F n | 4F 13 8F 13 | 8F 17 1QF 17 

s 2 s s s 2 s s 2 s s 2 s 

In these equations L = Infi/rrib, while the quantities F 1 ,...,F 2 i denote the following 
harmonic polylogarithms: 

F = [HPL({-2}, s), HPL({-1}, 1 - s), HPL({-1}, s), HPL({1}, s), HPL({2}, s), 
HPL({3}, s), HPL({4}, s), HPL({-2, 2}, s), HPL({-1, 2}, s), HPL({1, 1}, s), 
HPL({1, 2}, s), HPL({1, 3}, s), HPL({2, 1}, s), HPL({2, 2}, s), HPL({3, 1}, s), 
HPL({-1, 0, 0}, 1-s), HPL({1, 1, 1}, s), HPL({1, 1, 2}, s), HPL({1, 2, 1}, s), 
HPL({2,1,1},S),HPL({1, 1,1,1},$)] . (17) 

We now turn to the order a 2 contributions F)f\ which we decompose according to 



D?> = C F 



p(2) d(2) d(2) p(2) 

-"(-4),* , -"(-3),i , fVgjj -"(-!),» „(2) 

1 5 1 n 1 1" -TV, 



c4 ^ .3 ^ „2 ^ , ^-"-(0),i 



(18) 



The (infrared) singular pieces yield relatively compact expressions. We find 

d(2) ... 

K (-4),i - 2 
( 2 ) _ ri 5 xorxo^ , 11Ca 



R\% A = C F ^- + 2L + 2F 4 j+^-n l T R 



R?L 1 = C F — + 10L + 4L 2 + — + 8F 4 + 8LF 4 - -4 + 2F 5 + 8F 10 ) + 



,_ 2) ,, - „, v £ + W l + 4l» + £ « 

,'49 11L 7t 2 11F 4 \ 8 r _ / 10 41 4F 4 \ 
„f2) ^ /213 19C(3) 73L r2 16L 3 llvr 2 Ltt 2 55F 4 

leL *F t + ^p- - ^ - ^ + 8F 5 + SLF, - § + 2F 6 + 28F 10+ 
3 2s s s 



32LF 10 - ^ 



\ „ ( 1549 11C(3) 67L 
+ 8F n + 12F 13 + 32F 17 j + ^- — + -f> - 



7vr 2 Ltt 2 67F 4 tt 2 F 4 \ /20L vr 2 16LF 4 

+ ^ — 9- + ^J + (^ +8L + 

125 20L vr 2 20F 4 . 



p (2) _ R (2) _ n 

^(-4),2 ~~ (~3),2 ~~ U 

d(2) r f 2 , 2F 4 2F 4 



■<(-2),2 - _§ T- ~ 2 



(2) / 9 8L 7F 4 8LF 4 11F 4 8LF 4 2F 5 2F 5 12F 10 12F 10 

^-112 ~~ F I ~ — + — ^- H ^ z r — ^ 1 ^ — 

1 ; ' V s s s z s z s s s z s s z s 



L (-4),3 - (~3),3 



0(2) r I 2 2F * + AF * 



d(2 ) n /9 8L 7F 4 8LF 4 24F 4 16LF 4 2F 5 4F 5 12F 10 24F i0 

-Rf_il s — ~ H — : ^ ^ 1 : 1 : ttt H : tt, 1 : — 

v ; ' \s s s s z s s s s s z s 

(2) 

On the other hand, the expressions for the infrared finite parts i are rather lengthy. 
It is convenient to further decompose them according to 

„(2) \- °Fft,j,k + c Afj,j, k + n i T nfS,k + n hT R fiXk 
n (o),i 2^ §J(l- s) k 

In the following we list the functions fc , fi,j,ki fij,k an d fij,ki (2 = 1, 2, 3) for all values 
j, k for which they are nonzero. We find 

1327 16C(3) 213L 76C(3)L „ or2 80L 3 16L 4 97tt 2 „ 2 

f? n n = + V ; + + 73L 2 + + + 41n 2 vr 2 + 

Jl '°'° 16 3 2 3 3 3 48 V ; 

llLvr 2 2L 2 tt 2 449tt 4 4tt 2 F 1 10tt 2 F 3 153F 4 28C(3)F 4 „_ rrn 

1 H H sv ' + 110LF 4 + 

3 3 720 3 3 2 3 

2 64L 3 F 4 10tt 2 F 4 4 r 2 19F 5 2 tt 2 F 5 
64L F 4 H 1 V -Ln F 4 — + 32LF 5 + 16L F 5 + — 

12F 6 + 8LF 6 - 6F 7 - 16F 8 + 40F 9 + 59F 10 + 112LF 10 + 64L 2 F 10 + — — ^ + 

o 

28F n + 32LF n - 8F 12 + 60F 13 + 48LF 13 + 12F 14 + 12F 15 + 104F 17 + 128LF 17 + 
32F 18 + 48F 19 + 56F 20 + 128F 21 

fllfi = ^-^- 26LFl - iL i Fi + ^-^_ iLFs + Ft + m _ 

6 2 o I 

25F 10 - 24LFi - AF n - 10F 13 - 28F 17 



f 1,0,1 = -30C(3) + ^- + 161n(2)7r 2 + — + vr 2 F 2 ^ + + 90F 5 - 

4tt 2 F 5 + 12F 6 + 8F 7 - 80F 9 + 50Fi + 2AF U - 78F 13 - 8F 14 + 16Fi 5 + 2F 16 

7i,o,2 3 90 3 3 3 5^ 3 ^6 

20F 7 - 32F 8 + 32F 9 - 56F n + 112F 13 + 52F M - 104Fi 5 

fi,o,s = ^ + ^-3^-^ + 6^ + 24^ + 32^-6^-56^ + 
112Fi5-6Fi 6 
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89437 19C(3) 3925L 299L 2 44L 3 815tt 2 

-T296" + "18 54~ + 22C(3)L " -9 9 216" + 

16Ltt 2 2L 2 tt 2 31tt 4 2tt 2 F 1 5vr 2 F 3 2545F 4 A/ 466LF 4 

+ + + + 14C(3 )F 4 - 

9 3 120 3 3 54 sv ; 9 

44L 2 F 4 28tt 2 F 4 4 r „ 116F 5 44LF 5 4vr 2 F 5 20F 6 

„ 349F 10 88LF 10 4tt 2 F 10 UF u „ 62F 13 88F 17 

20F 9 + + 8Fi 2 

9 3 3 3 3 3 

1 2rn 269F 4 22LF 4 2vr 2 F 4 11 F 5 13F 10 

-r Fa+ V + ^-V + V- 4F9+ ^ 

15 C (3) + 13^ - W + 3 4-?Il + ^ + ^ + nFs- 



5 2 3 6 

' ' ' r i6 



4tt 2 F 5 - 31F 6 + 8F 7 + 40F 9 + 17F 10 + 13F n - 11F 13 - 8F 14 + 16F 15 - F : 



67tt 2 86vr 4 4tt 2 F 1 4vr 2 F 3 29vr 2 F 4 _ 38vr 2 F 5 

+ -^r 2 ^ 1 - 17i? s + — =— ^ + 30F 6 



6 45 3 3 3 

36F 7 + 16F 8 - 16F 9 - 14F n + 28F 13 + 20F 14 - 40F 15 

263tt 4 4tt 2 F 1 3n 2 F 2 29vr 2 F 5 ^ ^ ^ 

3-^ + — 3-^ - 3F 6 + 30F 7 - 16F 8 + 3F 13 - 14F 14 + 

28F 15 + 3F 16 



6629 26C 3 682L 100L 2 16L 3 85tt 2 8Ltt 2 418F 4 

I H 1 1 1 1 1 - + 

324 9 27 9 9 54 9 27 

152LF 4 16L 2 F 4 8tt 2 F 4 76F 5 16LF 5 8F 6 152F 10 32LF] 

9 + 3 + 9 + "9" + 3 + ~3~ + 9 + 3~ 
16F n I6F13 32F 17 

3 3 3 

38F 4 8LF 4 4F 5 8F 10 

9 3 3 3~ 



7951 28C(3) r r9 112L 3 41vr 2 2Lvr 2 530F 4 

s v ' +16L + 20L 2 + — — + — — + — — - + 8LF 4 + 



162 9 9 54 3 27 

r2r , 2vr 2 F 4 76F 5 I6LF5 8F 6 
16L 2 F 4 + ^ + + + 



9 9 
38F 4 8LF 4 4F 5 



9 3 3 

508 64tt 2 440F 4 104F 5 



+ 



^ 9 +16C( 3) + ^ + ^_48F 5 -16^ 
-16C(3) - + ^ + 18F. 
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-31 - 36L - 16L 2 + 3tt 2 + - 24F 4 - 44LF 4 - 16L 2 F 4 - 19 ^ ^ 4 + 

23F 5 - 8LF 5 + 18F 6 + I6F9 + 50F 10 - 48LF 10 - 24F n + 12F 13 - 56F 17 

-^tt 2 F 3 + I6F4 + 28LF 4 + 16L 2 F 4 - + 13F 5 + 8LF 5 - 2F 6 - 16F 9 - 

2F 10 + 48LF 10 + 8F n + 20F 13 + 56F 17 
8F 10 

_ 4C (3) + ?^ - 2,*F 2 + ^ - ™p + 44F 5 + 32F 9 + 100F 10 - 

24F n + 44F 13 - 4F 16 

w . 118tt 2 12tt 4 „ 16vr 2 F 3 136vr 2 F 4 „ „ 

32C(3) - — — - 8vr 2 F 2 + — ^ ^— * - 100F 5 + 16tt 2 F 5 + 

32F 6 - 32F 7 + 64F 9 - 112F n + 208Fi 3 + 32F M - 64F 15 - 16Fi 6 

304tt 4 16n 2 F 1 „ 136tt 2 F 5 ^ A „ 

—— + — 5-^ - 6tt 2 F 2 — — + 12F 6 + 48F 7 + 64F 8 - 12F 13 - 

45 6 6 

112Fi4 + 224Fi 5 - 12Fi 6 

269 44L n 2 2tt 2 F 3 257F 4 44LF 4 10vr 2 F 4 46F 5 A ^ o ^ 

+ 2tt 2 + + - + + 4F 6 - 8F 9 + 

9 3 3 9 3 3 3 

86F 10 

2tt 2 F 3 203F 4 44LF 4 4vr 2 F 4 22F 5 26F 10 
~3 ~9 3~ 2 + ^^^ + 8F9 "^ 

2C(3) + ^ + tt 2 F 2 - - + 32F 5 + 10F 6 - 16F 9 + 34F 10 + 

2F n — 2F 13 + 2F 16 

, 67tt 2 6tt 4 a „ 8vr 2 F 3 58tt 2 F 4 „ 
-16C(3) - — + 4vr 2 F 2 — ^ — i - 34F 5 + 8vr 2 F 5 + 

68F 6 - 16F 7 - 32F 9 - 28F n + 64F 13 + 16F 14 - 32Fi 5 + 8Fi 6 
263tt 4 87r 2 F! 9 58tt 2 F 5 

— + 3tt 2 F 2 — ^ - 6F 6 + 60F 7 - 32F 8 + 6F 13 - 28F 14 + 

yu o o 

56F 15 + 6F 16 

1 76 16L 52F 4 16LF 4 _ 8F 5 16F 10 

•/2,i,o - -y - -g g 3 3 3— 

, 52F 4 16LF 4 8F 5 16F 10 

/2,2,o - 9+3+3+3 



76 16L 292F 4 16LF 4 8F 5 



9 3 9 



rnh 
■'2,1,0 

,„„ 52F 4 16LF 4 s/;-, 

72,2,0 — n o 



9 3 3 



h _104 32tt 2 80F 4 16F 5 

/2,o,i 3 9 3 3 

, 64tt 2 32F 5 

72,0,2 - g 3 

/ 2 n o,3 = -32C(3) + 32F 6 
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31 + 36L + 1QL 2 - 3tt 2 - + 75F 4 + 9QLF 4 + 32L 2 F 4 + HZLZi _ 

12F 5 + 16LF 5 - 20F 6 - 32F 9 + 96LF 10 + 32 F n + 8Fi 3 + 112F 17 

-^-^ - 16F 4 - 28LF 4 - 16L 2 F 4 + - 13F 5 - 8LF 5 + 2F 6 + 16F 9 + 

2Fi - 48LFi - 8F n - 207\ 3 - 56Fi 7 
— 8F 10 

-40 C (3) - ^ + 161n(2 y + 4 ^F 2 - ^ - 4F 4 + ^ + 12F 5 - 

8F 6 - 64F 9 - 48Fi + 8F n - 56Fi 3 + 8F i6 

>x . 248tt 2 12tt 4 - 8tt 2 F 3 1167r 2 F 4 „ 
12C 3 - — — + 6vr 2 F 2 + — - -° — -i - 284F 5 + 16tt 2 F 5 + 

32F 6 - 32F 7 + 32F 9 - 300Fi - 88F n + 188F 13 + 32F 14 - 64F i5 + 12F 16 

-64C(3) + IIStt 2 + + + 16vr 2 F 2 - 16n 2 F 3 + 136tt 2 F 4 + 

45 3 

3OOF5 - 32 °^ ^ 5 - 112F 6 + 128F 7 + 128F 8 - 192F 9 + 336F n - 
640Fi 3 - 256F U + 512Fi 5 + 32F 16 

_304^_ _ lg7r 2 Fi + 187r 2 F2 + 1367r 2 F5 _ 36F6 _ 144 p 7 _ ig2Fg + 36Fi3 + 
15 

336F U - 672Fi 5 + 36F 16 

269 44L n 2 4vr 2 F 3 592F 4 88LF 4 14tt 2 F 4 68F 5 „ „ 

■ + 2tt 2 + — — i — i + — i - 4F 6 + 



9 3 3 
16F 9 - + 4F n - 8F 13 



2 „ 203F 4 44LF 4 4vr 2 F 4 22F 5 
--tt 2 F 3 + — — - + — — - — i + — - - 8F 9 + 



26F 



10 



9 3 3 3 3 



-6 + 20C(3) + 8\n(2)n 2 - 2n 2 F 2 + ^ - 18F 4 + ^± - 16F 6 + 

32F 9 - 4F 10 + 16Fn - 12Fi 3 - 4F 16 

A/ , 118tt 2 12vr 4 „ 4tt 2 F 3 89tt 2 F 4 1 1 ^ „ 
-6C(3) — - 3tt 2 F 2 — ^ — ^ - 144F 5 + 16tt 2 F 5 + 



3 5 3 3 ' 1 

82F 6 - 32F 7 - 16F 9 - 102F 10 - 54F n + 102Fi 3 + 32F U - 64F 15 - 6F 1( 

9 398tt 4 167r 2 F! , , 9 176tt 2 F 5 

32C(3) + 67tt 2 + — - 8tt 2 F 2 + 8tt 2 F 3 + 58tt 2 F 4 + 102F 5 — 

45 3 3 

196F 6 + 160F 7 - 64F 8 + 96F 9 + 84Fn - 184F 13 - 96F M + 192Fi 5 - 16Fi 6 
2637T 4 

+ 87r 2 F! - 9tt 2 F 2 + 58tt 2 F 5 + 18F 6 - 180F 7 + 96F 8 - 18F 13 + 84F 14 - 



30 

168F 15 - 18Fi 6 



nl 76 16L 152F 4 32LF 4 16F 5 32F 10 

/3,i,o - "g" + ^ _ + — ~ + 3 + ^"~ 

nl J2F\ _ 16LF 4 _ 8F5 _ I6F10 ' 

7 3 ,2,o - 9 3 3 3 
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h 76 16L 392F 4 32LF 4 16F 5 

, ' 3 ' 1 ' ~ "9 _ + ^ _ + — 9" + 3 + 

h 52F 4 16LF 4 8F 5 

h,2,o 9 3 3 

h 32 16tt 2 80F 4 32F 5 

/3,o,i ~~ g 9 3 3 



3 3 3 

320vr 2 352F, 

y 



/s,o,3 = - 64 C(3) — + —5— + 64F 6 



fsAA = 96C(3)-96F 6 . 



4 Two-loop results for the Wilson coefficients d 



In this section we give results for the Wilson coefficients Cj, valid to NNLO in a s . To 
calculate them, we take the UV-renormalized form factors Di obtained in the previous 
section, translate them to four-flavor QCD using ( Tl3|) . and evaluate the matching condi- 
tion ( flOi) . This procedure allows us to determine both the Wilson coefficients G{ and the 
renormalization factor Zj. The form of the renormalization factor is completely deter- 
mined by the renormalization-group equations for heavy-to-light currents in SCET, and 
thus provides important checks on our result. We shall first say a few words about these, 
and then list results for the Wilson coefficients Cj. 

The renormalization factor Zj is determined from our calculation by requiring that 
the matching relation ( flOl) is finite in the limit e — > 0. However, as explained in Section 
|2l it can also be determined by the UV poles of the object J ® S in the parton model. 
Expressions at two loops can be derived from the renormalization factors for the jet and 
soft functions calculated in [14,16]. Either way, the result depends only on logs of the 
form L p = In /i/n • p and reads 



Air 



1 1 

"- + - 

e z e 



5 
2 



2L, 



Air 



r(2) 



(19) 



where the two-loop coefficients are 



7(2) 
'(-4) 

7(2) 
V (-3) 

7(2) 
\-2) 



7 (2) 
Z (-l) 



2 

25 



GVI — + r >L 



2L: 



niT R 



3 n z 



+ CU(— + — 



M 54 6 9 



18 ' 12 
1549 



216 



11L, 



III 
24 



+ n z T fl - — 



11 



10 


4L p 




3 




7T 2 " 


y 4 


~ y 



(20) 



In SCET, the hard function is derived from the matrix of Wilson coefficients = CiC. 
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and satisfies the renormalization-group equation [4] 
d 



dhifj, 



Hij(n-p,ii) = 2 



i(a s ) + r cusp (a s )ln 



n ■ p 
H J 



Hij(n-p,fji) . 



(21) 



It is easy to show that the anomalous dimension derived from the explicit expressions in 
201) are consistent with (ED), with 

1549 7tt 2 11 

A 1 + — r C3 



7 



4tt V47t/ [ V 16 4 7 



432 48 



,125 7T 2 
-TIiTr 1 

1 108 12 



'.22) 



This result is consistent with that given in [33], and the piece of the anomalous dimension 
proportional to the logarithmic term is consistent with the two-loop cusp anomalous 
dimension from [34]. 

We now give final results for the Wilson coefficients Ci, which we decompose according 

to 

a = cf + ^cf + (^)V>. P3) 



We find 



C[ 0) 



1, cf ] = cf ] 







(24) 



a 



(i) 



R [ V. , Cf } = R™. + C 2 F K hl + C F n h T R K it2 (i = 1, 2) , 
1- s 



-R 



(i) 



a 



(2) 



1 - S 



R 



(2) 



+ ClK^ + C F n h T R K 3 o 



2 (0),3' ^3 2 V'"'^' 3 

where the R$,j were given in Sectional and functions and K i:2 read 
49tt 



-54- 
13F 5 



^ 7T 4 



60F 4 



13tt 2 F 4 I8F4 7T 2 F 4 31F 5 7T 2 F 5 
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+ 



+ 



12s 



48 160 

^ - 8F 6 + -4 - 2F 7 - 63F 10 - tt 2 F 10 + - 22F n + 

2s s s s 



2 6 



+ 



6 P 1 fiF 

8F 12 - 28F 13 + — - - 8F 14 - 12F 15 - 68F 17 + — — - 24F 18 - 32F 19 - 

s s 

32F 20 - 80F 21 + ^ + ^F 4 C(3) + 



-7S 



5tt^ 



71F 4 - vr 2 F 4 + 



21F 4 4F 5 16F 10 

— — - 22F 5 + - 8F 6 - 68F 10 + — — - 

s s s 



24Fii - 32F 13 - 80F 17 + 



4C(3) \ 



97 7T 2 

( -— - — -38F 4 + 
1 2 2 s 



3 J 
6F 4 



12F 5 -40F 10 \+L i - — 



50 40F 4 \ 10L 



5tt 5 



2tt 2 F, 



L 2 (20 + I6F4) - 



9 

112L 3 
9 



-16 



2vr 2 



-8F 4 + 



8F 4 16F R 32F, 



10 



3s 



14 



18 7T 2 9F 4 7r 2 F 4 17F 4 tt 2 F 4 7F 5 7F 5 2F 6 2F 6 22F 10 

s 6s s z bs z s bs s z s s z s s z 
22F in 8F11 8F11 12F 13 12F 13 32F 17 32Fi 



+ — ; ^ + — — 1 ^~ + — + 



s s 2 s s 2 s s 2 s 



26 18F 4 + 26F 4 _ 8F 5 + 8F 5 _ 32F 10 + 32Fi \ + L 2 f 12 _ 12F * + 12F * 



s 2 s s 2 s s 2 s J V s s 2 



16 16F 4 16F 4 

A22 = L 7. — 

1 3s 3s 2 3s 



18 7T 2 9F 4 tt 2 F 4 53F 4 tt 2 F 4 7F 5 20F 5 2F 6 4F 6 22F 10 

3, 1 n ^ ~ ' ^9 : o 9 -"•9 A 9 

s 6s s^ 6s z S OS S 8 S 8 s z 

80F 10 8Fn 16Fn 12Fi 3 24F 13 32F i7 64iq 7 



s' 2 s s 2 s s 2 



S" S S" S S" s 

26 18F 4 68F 4 8F 5 16F 5 32Fi 64Fj 

s s s s z s s s 

' 12 12F 4 24F 4 \ 



s s 2 



r 16 16F 4 32F 4 

A 3 9 = L 

1 3s 3s 2 3s 



The terms proportional to the explicit factors of rih in (12411 stem from converting the 
results of the renormalized form factors Di from the five-flavor to the four-flavor theory. 
As a final check, we have confirmed that the /i-dependence in the Cj is such that the 
renormalizat ion-group equation (|2ip is satisfied. 



5 Conclusions 

We have presented results for the short- distance Wilson coefficients needed to complete 
the calculation of partial decay rates in B — > X u £u at NNLO in a s and to leading order in 
1/mb, for decay kinematics limited to the shape-function region. The technical challenge 
was to compute the two-loop QCD corrections to the semi-leptonic b — > u transition 
current. To do this, we used the Laporta algorithm to perform a reduction to master 
integrals, which were solved using the method of differential equations. We then performed 
a matching calculation from QCD onto SCET to translate these results into the Wilson 
coefficients needed to compute the hard function in the factorization formula (CEJ) at NNLO. 
In a companion paper, we shall perform an analysis of partial decay rates with arbitrary 
kinematic cuts at NNLO, and study the implications on the determination of \V u b\ from 
inclusive decays. 
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Note Added: After our calculation was completed, the paper [35] appeared, where 
the UV-renormalized two-loop corrections to the b — > u current were presented. We have 
compared with their results and found agreement with those given in Section [3l 
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